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In this paper the effect of a small dissipation on waves is included to find exact solutions to the 
modified BBM equation. Using Lyapunov functions and dynamical systems theory, we prove that 
when viscosity is added to the BBM equation, in certain regions there still exist bounded traveling 
wave solutions in the form of solitary waves, periodic, and elliptic functions. By using the canonical 
form of Abel equation, the polynomial Appell invariant make the equation integrable in terms of 
Weierstrass p functions. We will use a general formalism based on Ince's transformation to write the 
general solution of dissipative BBM in terms of p functions, from which all the other known solutions 
can be obtained via simplifying assumptions. Using ODE analysis we show that the traveling wave 
speed is a bifurcation parameter that makes transition between different classes of waves. 

PACS numbers: 02.30.Gp, 02.30.Hq, 02.30.Jr 



I. INTRODUCTION 

In 1895, Korteweg and de Vries, under the assump- 
tion of small wave amplitude and large wavelength of in- 
viscid and incompressible fluids, derived an equation for 
the water waves, now known as the KdV equation [9], 
which also serves as a justifiable model for long waves in 
a wide class of nonlinear dispersive systems. KdV has 
been also used to account adequately for observable phe- 
nomena such as the interaction of solitary waves and dis- 
sipationless undular shocks. For the water wave problem 
(LI) is nondimensionalized since the physical parameters 
\u\ = 3r]/2H, X = V6x7i^, and t = ^J^g/Ht*, where r] 
is the vertical displacement of the free surface, H is the 
depth of the undisturbed water, is dimensional dis- 
tance and the dimensional time are all scaled into the 
definition of nondimensional space x, time and water 
velocity u{x,t). When the physical parameters and scal- 
ing factors are appropriately absorbed into the definitions 
of ix, X and t, the KdV equation is obtained in the tidy 
form 



- UUr 



0. 



(LI) 



Although (LI) shows remarkable properties it mani- 



fests non-physical properties; the most noticeable being 
unbounded dispersion relation. It is helpful to recognize 
that the main difficulties presented by (I.l) arise from 



the dispersion term and arise in the linearized form 



Ut 



0. 



First, note that when the solution of (1.2) is express- 



ible as a summation of Fourier components in the form 
F(/c)e~**^^^+^^\ the dispersion relation is ijj(k) = k — k^. 
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The phase velocity Cp = becomes negative for 

/c^ > 1, contrary to the original assumption of the for- 
ward traveling waves. More significantly, the group ve- 
locity Cg = ^ = 1 — 3k^ has no lower bound. To cir- 
cumvent this feature, it has been shown by P eregrine [13], 
and later by Benjamin et. al. [3 that ( |I.1[ ) has an alter- 
native format, which is called the regularized long- wave 
or BBM equation 



Uxxt 







(1.3) 



where the dispersion term Uxxx is replaced by —Uxxt to 
refle ct a bounded dispersion relation. It is contended that 
is in important respects the preferable model over 
which is unsuitably posed model for long waves, 
he initial value problem for this equation has been 
investigated previously by [TSHIS]. In [3], the authors 
showed that its solutions have better smoothness prop- 
erties than those of (I.l). BBM has the same formal 



justification and possesses simil ar p roperties as that of 
KdV. The linearized version of (1.3) has the dispersion 
relation ijj(k) = j^^, according to which both Cp, and 



Cg are bounded for all k. Moreover, both velocities tend 
to zero for large /c, which implies that fine scale features 
of the solution tend not to propagate. The preference 
of (1.3) over (I.l) became clear in 13 , when the authors 



attempted to formulate an existence theory for (LI), re- 



spective to nonperiodic initial condition u{x,0) defined 
on (— oo, oo). 



(1.2) provided by the equation 



A generalization to (1.3) to include a viscous term is 



Uxxt 



(1.4) 



where > is transformed kinematic viscosity coefficient 
of a liquid is a popular alternative to the Korteweg-de 
Vries equation, for modeling unidirectional propagation 
of long waves in systems that manifest nonlinear and dis- 
persive effects. In [3 , the equation is recast as an integral 
equation, and solutions, global in time, are found for ini- 
tial data in ly^'^(R). In [5 , global solutions are found 
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for initial data in i7^(R), 5 > 0. In this article, we study 
the effect of dissipation on (1.3), and henceforth we shall 



refer to (1.4) as the modified BBM by viscosity. 



II. EXISTENCE OF DISSIPATIVE SOLUTIONS 

The existence and stability of solitary waves has been 
investigated previously by Bona et. a/., lElHll^. In this 
section we propose an easy theorem based on dynamical 
system theory, in which we prove the existence of solu- 
tions to (1.4) for both c < 0, and c > traveling wave 



velocities. Also, we will show that the solution tends to 
the fixed points of the dynamical system c — 1 ± /c, for 
any k and large (. So, at the both ends of the boundaries 
the solution does not necessarily decay to zero, but there 
is fiow of energy through the boundaries to the solution 
and back to the boundaries. 

Theorem 1 Let c 7^ and k > 0. ( |I.4| ) has a traveling 
wave solution of the form u{x^ t) = u{x — ct) = u{Q with 
u{() {c — 1) -\- k as ( ^ —00 and u{Q {c — 1) — k 
as ( ^ 00. If c < 0, then u{Q < {c — 1) -\- k for all real 
(. If 0, then u{Q > {c — 1) — k for all real (. 

Proof 

The class of soliton solutions is found b y em ploying the 
form of the traveling wave solutions of (1.4) which takes 
the form of the ansatz 



u{x,t) = u{C) 



(II.l) 



where ( = x — ct is the traveling wave variable, and c is 
a n onze ro translational wave velocity. The substitution 
of (|TlT| in (O) leads to 



(1 - c)ii' + W + cu''' = uu'\ 



(II.2) 



where ' = 4?^. Integrating, 



(1 - c)u + -u^ + cu" = uu' (II.3) 



for some constant A. With cj) = u and ifj = u' ^ (II.3) is 
equivalent to the dynamical system 

^' = ^ 

c c Zc c 



k —(c—1) 

for A = ^ — —. We will show that (II.4) has a trajec- 
tory from z+ to z~. 

First, we will consider the c < case. The linearization 



of (II.4) at z+ is 



1 " 




> - 0+ " 




> - 0+ " 


^(c-l)2+2A ^ 

c c _ 






= M 





(11.7) 

The determinant of M is negative, so z+ is a saddle point, 
and M has eigenvalues with A~ < < A+. 

The unit eigenvector corresponding to A+ and pointing 
to the left is 



-1 
-A+ 



yiTTA+p' 



(11.^ 



Letg(C) 
z+ and 



51(c) 
52(c) 



be a solution of (II.4) with g(C) 



g(C) 



|g(C) 



(11.9) 



as C ^ —00. We will show that g(C) ^ z as C ^ 00 and 
9i{C) < foi" all I'^al (. Define 



P3(<^) = ^<^' + ^(l-c)(/.2-A,^. 



(11.10) 



P3 is increasing on (— 00, (/) ], decreasing on[<p ^4''^]^ and 
increasing on [^+,00). Define 



£(u) =Pj(<A)-P3(0+) 



1 



!.2-^0-P3(,/.+ ) + -C^2_ 

(11.11) 



Then for all C, 

|£(g(C)) = (Ig! + (1 - c)g^ - A)g[ + cg^g'^ = 
= (^5i + (1 - c)gi - A^g2+ 

■■i— 

\ c 



1 



V A\ 2 

-^2 + - = > 0- 

c c J 

(11.12) 

Since /:(z+) = then C{^{0) > for ah real (- 
Let 6 < 6~ with 



") < ^^3(0+), 



(11.13) 



u' = F(u) 



in vector notation. The fixed points of (II. 4) are 













(11.5) 



c-l± ^(0-1)2 + 2^ 




' c-l±k' 














(II.6) 



and let 

B = Ps{cj>-)-Ps{<P+) = max{Ps{<j>) U~ < < <P+}-P3{<j>+). 

(11.14) 

Let L > be large enough so that 



1 



cL^ < 0, 



(11.15) 
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and define 



(11.16) 



We claim that g(C) G R for all real As C ^ — 
g(C) z+ along the direction [1,A+]^. Therefore 
>^(g(C)) > ^(2:^) = for all real ("^ and there exists 
C G M such that g(C) G for ah C < C- 

So, it suffices to show that £ < on dR. Let u G 
If = then i2(u) = ^c?/^^ < 0. Similarly, if = 
then C{m) = P^{(^—) - Psict^^) + ^c^^^ < 0. If (/)— < 
^ < ^+ and = L, then jC{u) = P^{(^)-P^{(^+)^\cL^ < 
B + \cL'^ < 0. Similarly, if (j)— <</)<(/)+ and = -L, 
£(u) < 0. 

Define 



hm C{g{C)). 



(11.17) 



I/oo exists and is finite because £(g(C)) is a non- 
decreasing function of is bounded on R, and g(C) G 
R for all real Let C be the cj-limit set of g, i.e.. 



C = {u I 3(Cm) C M with lim 



00, g(C„ 



^u} 
(11.18) 

C is nonempty because g(C) G for all real ( and is 
bounded. C is compact and connected, and C = L^o on 
C. We must prove C = {z~}. 

Suppose (7 contains a point w that is not on the (/)-axis. 
Then w is not a stationary point of (II.4), and C contains 



another point y. Let ^ > with S < min(|w — y|, \w2\)/2. 
Then g crosses the annulus ^^(w) \ 5^/2 (w) infinitely 
many times, and there exist sequences (Cm) 7 (Cm) with 
Cm ^ 00 as m ^ 00, Cm < Cm < Cm+1 for all m, 
|g(Cm) - w| = (5/2, |g(C'm) - w| = (5 and S/2 < |g(C) - 
w| < (5 for all m and all C G (Cm, Cm)- 

Let = sup(|F(u)| | u G ^^(w)) < 00, where F is 
from dlLSf . Then Cm - Cm > S/{2K) for ah m. For 

C ^ (Cm 7 Cm)? 



^/:(g(C))>.^i(C)>^>o, 



(IL19) 



> 



But Loo is finite. This is a contradiction. 

Therefore, C is a subset of the 0-axis. Since C is 
connected and compact, C = [^1,^2] x {0} for some 
< 4>i ^ 4>2 < (t^'^ ' Recall that C = Loo on C. C 
is not constant on any segment of the ^-axis of positive 
len gth. So C is a single point, which is a stationary point 
of ( |II.4[ ) in R. There is only one such point, z~. 

Now consider the case c > 0. Let /c > 0. Since — c < 0, 



we have seen that there exists g solving (II. 2) with c 
replaced by — c, that is. 



(l^c)g' ^gg' -eg' 



^9 



0, 



(11.22) 



with g{C) (-c-l) + /c as C ^ -00, g{() {-c-l)-k 
as C ^ 00, and ^(C) < (— c — 1) -\- k for all real C- 

Define : R ^ R by u{C) = -g{-C) - 2. Then 
9(0 = -u{-C) - 2, ^'(C) = u\-C). 9"{0 = -u"{-0. 
and g"'{C) = "^'"{-Q- Substituting into ( |II.22| ) yields 

= (1 + c)u\-0 + {-u{-0 - 2)u\-0 - cu"\-0 + vu'\-0 

= (c - l)u\-C) - u{-C)u\-C) - cu"\-C) + vu'\-Q 

= -[(1 - c)u'{-0 + u{-Ou'{-0 + cu"'{-0 - vu"{-0]. 

(11.23) 

Therefore u satisfies (|II.2|). By definition of 



lim u{0 = lim (-5(-C) - 2) = lim (-5(C) - 2) 

C^ — cx) C^~c>o 



-2 - lim 5(C) 
(c- l) + fc, 



and 



-2-((-c-l)-fc) 

(11.24) 



lim u{Q 



lim (-^(-0- 2)= lim (-5(0-2) 



-2- lim 5(C) 

C^— 00 



= (c-l)-fc. 



((-c-l)+fc) 



(11.25) 



For all real C, g{C) < (— c — 1) + fc, so for all real C, 

^^(C) = -5(-C) - 2 > -[(-c - 1) + fc] - 2 = (c - 1) - fc. 

(11.26) 

Theorem [1] is proven. 



^{Cm - Cm)^2 ^ 



JII. EXACT SOLUTIONS VIA p FUNCTIONS 



> 



8K 



for all m. Therefore 

Loo=>C(g(Cl)) + (ioo-Ag(Cl))> 

00 

> £(g(Cl)) + ^ C{g{Cm)) - C{g{U)) > 



(11.20) 



A. Abel's equation 



(11.21) 



> 



Ag(Ci)) + E 



m=l 



= 00. 



The importance of Abel's equation in its canonical 
forms stems from the fact that its integrability leads to 
closed form solutions to a general nonlinear ODE of the 
form 

^CC + Qo{u)u(^ + q2{u) = . (ULl) 
This can be expressed by the following lemma. 
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Lemma 1: Solutions to a general second or- 
der ODE of type (III.l) may be o btaine d via 
the solutions to Abel's equation ( IIL6| ), and 
vice versa using the following relationship 



(IIL2) 



Proof: To show the equivalence, one just need the 
simple chain rule 



dr] 



drj 



(III.3) 



which turns (|IIL1) into the second kind of Abel equation 
-^ + goM^ + g2M =0. (IIL4) 



drj 



du 

Moreover, via transformation of the dependent variable 

1 



r]{u) 



(IIL5) 



(|IIL4|) becomes the Abel first kind without free or linear 
dy 



terms 



du 



(IIL6) 



Comparing ( IIL1[ ) with ( |IL3[ ) we identify the constant, 
and quadratic polynomial coefficients 



Qo{u) 



1-c A 



(IIL7) 



Q2{u) = TT + ^ 

2c c c 



The need of dissipation is imperative since without the 



dissipation term qo = 0, the Abel'e equation (IIL6) be- 
comes separable, as we will see next. 



Progress of integration of (III. 6) is based on the linear 
transformation v = J q o{u)d u = —^u, which allows us to 



write Abel's equation ( |III.6 ) in canonical form 



(III. 



g{v) is the Appell's invariant and is the quadratic poly- 
nomial 



9{v) = 



c(l-c)^ 



A 



= a2'^^+ai'U+ao (III. 9) 



By Lemke's transformation [TT] 



y 



1 dz 

z dv 



(III.IO) 



(|III.8|) can be written as a second-order differential equa- 

riant g{v). 

+ g{v) = (III.ll) 



tion with a preserved invariant g{v). 



,d^v 



dz^ 



B. No viscosity 



If = (II.3) becomes 



The classical solutions obtained before by assuming 
A = are 



u{x,t) = 3(c- l)Sech^ 



rV(c-i)/c 



(x — ct) 



-Sec^ 



1 



2 



(x-t/(l + c)) 



if c> 1 
(III.13) 

if < c < 1 



see left and middle panels of FIG. 1. The solitary waves 
that move with velocity c > 1 will be called fast waves, 
while the periodic solutions that travel with velocity < 
c < 1 slow waves. 



In (|III.6|) if let = 0, then the Abel's equation 

Q2{u) 



dy 
du 



(III.14) 

is separable and leads to the elliptic algebraic equation 
1-c 



T] 



u 
3^ 



2A 

c 



2D, 



(III.15) 



which after the substitution u^ = r] becomes the ellip- 
tic differential equation with solutions in term s of th e 
Weierstrass p functions reducible to sech or sec (III. 13). 



Yc 



1 



2A 



2D. 



(III.16) 



For verification, if one takes ^ of the above, we recover 



(III. 12). We will not discuss the reduction to hyperbolic 



or trigonometric functions, since this was done extensible 
before by many authors, but instead we refer the reader 
to [12 . 

Instead, we will analyze the novel case of the solutions 
that travel with a critical velocity c = 1. At the bound- 
ary between the fast and slow waves we will encounter 
periodic solutions in terms of rational combination of Ja- 
cobian elliptic functions. 



If c = 1 and A = 0, eq. (III. 12) becomes 

,2 



U 

~2 



- u 



0. 



(III.17) 



By multiplying by u^ and integrating once we obtain 
where B is some nonzero constant of integration. In 



(III.18) 



(III. 15) if we let D 



we will also arrive to 
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(111.18). Now let's use the substitution (i^ = B — 



where fi = /i(C)- Hence (III. 18) is 



= - 3BAt^ + 3B2. (III.19) 



Moreover, let's assume 

m(C) = ^-/BziQ, then (|nLT9| 

becomes 



2^ 



1. 



(III.20) 



This differential equation will be solved using Jacobian 
elliptic functions. 



Claim 1 if Z = 1 cos 2a + z"^ , then 

r 

u{x^ k) = 

Jo 



dz 1 2x 
= -sn 



VZ 2 l+x2' 
with k = sin a, see Bowman 

Proof 

Putting z = tan6>, we find 



u{x^ k) 



dO 



/o \/l - sin^ a sin^ 26> ' 

followed by y = sin 2^, which will give us 

dy 



v/(i-y2)(i-fcV)' 

where the integrand is the elliptic integral of the first 
kind. Therefore, solving (III.20), we obtain the solution 
in z 



(III.21) 



2z , .VB . 



where k = sin y| = is the modulus of the Jacobian 

elliptic function. It follows that in /i we will have 

^^ = ±sn(^aC,k), (III.22) 

where a = \/B By solving the quadratic, we have 



l±cn(^aC,k) 
sn(^aC,k) 



(III.23) 



Therefore, the solution to the BBM eq. without vis- 
cosity and velocity c = 1 is 



=5 



sn(f(x-t),k) 
lTcn(f(x-t),k) 



(III.24) 



±cn(f(x-t),k)y 



These critical solutions are 2K periodic with K given by 
the complete elliptic integral 



K 

~2 



f 

Jo 



dz 



= 2.76806, 



(III.25) 



see right panel of FIG. 1. The same solution was previ- 
ously obtained by Cornejo-Perez and Rosu, using a fac- 
torization technique [16]. 

The expressions (|III.13|) and (|III.24|), describe the 



whole class of solitary wave solutions, with spectrum 
c G (0, oo) which is in concordance with [4]. They present 
periodic and solitary wave solutions to ( |I.3| ), with no re- 
strictions on the constant A, which may be advantageous 
when adapting results to physical problems. 

Remark 1: all the solutions obtained using simpli- 
fying assumptions v = , A = are particular cases of 
general solutions of (II.3) in terms of Weierstrass elliptic 
functions p, as we will see next. 



C. Viscosity present, jy > 



Consider ( |III.8 ) in the form of non- autonomous eq 
F{y^ yv,v) = 0. Painleve proved in 1887 that all integrals 
of non- autonomous equations do not have movable singu- 
lar points, but only poles and fixed algebraic singularities. 
Poincare [17 proved in 1885 that any non- autonomous 
equation of genus p = is reducible to Riccati eq, while 
if p = 1 is integrable via Weierstrass p functions, after 
a linear fractional transformation. Since Riccati equa- 
tion does not arise form a nonlinear ODE and the Appel 
invariant g{v) has no singularities. Therefore, once can 
conclude that the closed form solution will not have mov- 



able singular points, but probably poles in (III.ll ). Sinc e 
v{z) has poles of order two, then the solution to (III.ll), 
must be written in terms of Weierstrass p functions using 
the transformation 



EzPuj{z)^F, 



(III.26) 



see Ince [8 . By substituting this ansatz in (III.ll) with 



E, F constants, then if p = |, uj will satisfy the elliptic 
equation 



(III.27) 



with elliptic invariants ^2,^3- 

Moreover, z is an exponential function, because it sat- 
isfies 



which leads to 



dv 1 

dz yz 



c z 



(III.28) 



(III.29) 
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The p solutions can be combined into the general sub- 
stitution 



with invariants 5^2 = 0, and gs- If ^3 = 1, this is the 
equianharmonic case, see [1 . 



^0 = ^-^(0, 



(111.30) 



z{Q is found by integrating (III. 35) to get 



where g{() = e <j, and n are constants which 

are related to A a nd p. 

By substituting ( |III.3Q[ ) into (|lL3| we obtain the ODE 



Z{C) = C2 + 



5v3c uc 

p 5c 

6u 



(III.42) 



1-c + cr 



9 = 0. 



(III.31) 



The free term was eliminated by setting A = ^+(t(1— c). 
Since g{0 = e-"'^f}(C), ( |lIL3l| ) becomes 



nv 1 — c + (T^ 
c 



fi(C) = = p(c4 + ^e^^ , 0, 53), (III.43) 

where 04 = 02 + 03, and ^3 are two integration constants. 
Using all of the above, together with (III. 37) we obtain 



2c 



(III.32) 

Finally, let = and by chain rule with ' = ^ 

and ■ 



25c 



14^/2 //14i/2x2 



(III.44) 



^, (|nL32|is 



(zO'a)+(^'-(2n+^)z')cO+(n2 + 
By letting 



\ 1 

CJ = — t 

) 2c 
(III.33) 



Then, the general solution to (III. 27) is 



u(C) = 



/14l/2x2 



25c^VfeJ -2A-e-p(c4+ 

(III.45) 

see FIGS. 2,3. Finally, it is worth mentioning that the 
solution in terms of p functions is not unique. One 
can see that there is at least one other similar form 
previously obtained by Porubov, via the Carnevalle 
method and Backlund transformation [101 UHl • 



5V3c 



e^s0,^3), 



we obtain 



(III.34) 



(III.35) 



We also choose a = —{in? c + n^^ + 1 — c) which cancels 
the linear term in ( |III.33 ). 
We are left to solve 



2c 



Remark 2: L etting v = 0^n = 0, ci = l and 
c=l^cr = Oin ( III. 38 ), and sin ce uo = but / = — ^, 
then uj = —f so ( III.40p becomes (III. 17), with solution 
given by (|III.24[), see FIG. 1 right. 



(III.36) 



subject to (III. 35). Set n 



-1^, then 

5c ' 



Remark 3: Lett ing u = Q ^ n = 0, and a = c— 1 7^ 0, 
then according to (III. 35), z' = ci. Since uj = g^ then 
(III. 33) is actually (III. 12), and hence we recover the 



14^2 

1 . 

2 



slow and fast waves solutions (III. 13), see FIG. 1 left 
and middle. 



(III.37) 



A 



25c 



Remark 4: If the second invariant ^3 = in (III.40), 
then by integrating (III. 40) 



By substituting (III. 35) into (III.36), we obtain 



uj{z) 



2ccl 



(III.38) 



(III.46) 



Using (III. 42) in (III. 46) we obtain the kink solution 



Letting ci = we arrive at the elliptic equation 

ij = 6uj^ (III.39) 
which by multiplication by u and integration becomes 
cj^ =4uj^ -gs. (III.40) 



/(C) 



e 5c 



Its solution is 



uj{z) = p(2: + C3,0,^3) 



(III.41) 



(III.47) 

see FIG. 4. 

An analysis when a p function tends to a smooth kink 
can be found in [19 . 

Remark 5: lfA = 0^c-l = ±k = and one 
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selects the lower branch of the radical, we obtain 



/(C) = - e"^- p{c4 + ^\/^e^sO,53) if 53 7^ 



/(C) = 



2uC 

e 5c 



if gs = 0. 



(111.48) 



IV. STABILITY OF THE VISCOUS WAVES 

In this section we consider the two-mode dynamical 
system of ( |IL4[ ) with equihbrium points in the phase 
plane (0,?/^) at (0,0) and (2(c - 1),0), when A = 0. AW 
the equilibrium points lie only in the plane = 0, in 
the (^, i/j^ c) space. The bifurcation curves are given by 
0(0-2(c-l)) = 0, which are two constant lines that coin- 
cide at c = 1, which is a bifurcation point, since the num- 
ber of equilibrium points changes from two (0 < c < 1), 
to one (c = 1) and back to two (c > 1). Near the origin 
(0-,O) = (0,0) 



(IV.l) 



Following standard methods of phase-plane analysis the 
characteristic polynomial of the Jacobian matrix of 
(IV.l) evaluated at the fixed point ((/)~,0) is 



\4>'] 




" 1 ■ 




>" 






1 — c u 
c c . 







goiX) = -poX + qo = 0, 



(IV.2) 



1-c 



where Po = and qo 

Since > 0, c > 0, then po > 0, hence the origin is 

unstable . Also, putting A = p^ — 4go = ^, we 

have the following cases: 



(i) < c < 1, gives gi < ^ A > 0, hence the fixed 
point is a saddle point, see FIG. 5^ 

(ii) c > 1, gives qi > 0. if u > 2^ c(c — 1), we have 
an unstable node\ if < < 2^/c{c — 1), the fixed 
point is unstable spiral, see FIG. 6. 

All the other remaining cases, i.e, both fixed points 
collide, are degenerate, since if c = 1, ^(A) = — vX. 

Note that when there is no viscosity, the system has 
Hamiltonian, which is the Lyapunov function 



i*' = ic(u), (IV.5) 

be c 

see ( |II.11| ). Therefore, along any phase path H(u) = 
constant is a conserved quantity, see FIG. 8. In this 
case, po = pi =0, and hence 

(i) (^~,0) is a saddle and (^^,0) is a center when 
qi > 0, see FIG. 7 left, 

(ii) (0~,O) is a center and (0^,0) is a saddle when 
qi < 0, see FIG. 7 middle, and 

(iii) c = 1 is degenerate since both critical points collide 
at origin = see FIG. 7 right. 

Therefore, in this case the unstable spirals (which corre- 
spond to the case with small viscosity u = 0.1) become 
centers when there is no viscosity u = 0. 

This is an example of a transcritical bifurcation where, 
at the intersection of the two bifurcation curves ^ = and 
= 2(c — 1), the equilibrium changes from one curve to 
the other at the bifurcation point. As c increases through 
one, the saddle point collides with the unstable node, and 
then remains there whilst the unstable node or spiral 
moves away from (^~,0). 



V. SUMMARY 



(i) < c < 1, gives qo > 0. If u > 2y^c{l - c) > 0, the 
origin is unstable node; if < < 2y^c(1 — c), the 
origin is unstable spiral, see FIG. 5, 

(ii) c > 1, gives qo < ^ A > 0, hence the origin is a 
saddle point, see FIG. 6. 

Near the secondary fixed point (^+, 0) = (2(c — 1), 0) 



(IV.3) 



The c haracteristic polynomial of the Jacobian matrix of 
(IV.3) evaluated at the fixed point (^^,0) is 







■ 1 " 










1 — c u 
c c . 







^i(A) = A^-piA + gi =0, 



(IV.4) 



where Pi = Po = ^, and qi -qo = 

Since pi = Po-, the second fixed point is also unstable , 
and moreover we have the cases: 



In this paper a basic theory for the BBM equation 



(1.3), and its extension (1.4) to include the dissipation 



term was shown. When the viscous term is not present, 
(1.3) has traveling wave solutions that depend critically 



on the traveling wave velocity. When the velocity is sub- 
unitary, it has solitary wave like solutions. If the velocity 
is super unitary, the solutions become periodic. At the 
interface between the two cases, when c = 1, (1.3) has 



periodic solutions in terms of the Jacobi elliptic func- 
tions. Based on Ince's transformation, a theory based 
on the dissipative term was presented, in which we have 
found a set of solutions in terms of Weierstrass elliptic p 
functions. Based on dynamical systems theory we have 
proved that the solutions of ( |I.4[ ) experience a transcrit- 
ical bifurcation when c = 1, where at the intersection 
of the two bifurcation curves, stable equilibrium changes 
from one curve to the other at the bifurcation point. As 
the velocity changes, the saddle point collides with the 
node at the origin, and then remains there, while the 
stable node moves away from the origin. 
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Kink, c=1.5,v=l 




FIG. 4. Kink solutions zy = 1 eq. (III.47), left c = 1.5; middle c = 0.5; right c = 1 
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FIG. 5. Phase-plane for c = 0.5, left z/ = 1, (0.0) node, (-1,0) saddle; right v = 0.1, (0.0) spiral, (-1,0) saddle 

















































































































l-j U-Lu-,.- 














1 





























































































FIG. 6. Phase-plane for c = 1.5, left z/ = 1, (0.0) saddle, (1, 0) node; right v = 0.1, (0.0) saddle, (1, 0) spiral 
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Hamiltonkn, c=l .5 




Hamiltoiikn, c=0.5 



Haiiultoman, c=l 





FIG. 8. Contour plot of Hamiltonian H = const see (IV. 5), left c — 1.5; middle c = 0.5; right c—1 
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